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We describe a criterion for the detection of entanglement between two multi-boson systems. The 
criterion is based on calculating correlations of Gell-Mann matrices with a fixed boson number 
on each subsystem. This applies naturally to systems such as two entangled spinor Bose-Einstein 
condensates. We apply our criterion to several experimentally motivated examples, such as an S^S^ 
entangled BECs, ac Stark shift induced two-mode squeezed BECs, and photons under parametric 
down conversion. We find that entanglement can be detected for all parameter regions for the most 
general criterion. Alternative criteria based on a similar formalism are also discussed together with 
their merits. 
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I. INTRODUCTION 

The notion of quantum entanglement is closely con¬ 
nected to the tensor product structure of a state space 
of composite quantum systems. For some types of com¬ 
posite systems that lack this structure, entanglement is 
no longer uniquely defined. This happens for systems of 
indistinguishable particles, in which symmetrization or 
anti-symmetrization of the wave function imposes some 
sort of “default” entanglement. As an elementary exam¬ 
ple of this, consider two identical bosons which occupy 
two quantum states, for example atomic levels denoted a 
and b of an atom. As they are indistinguishable bosons, 
it is appropriate to use the second quantized notation to 
write the quantum state as 

a^6t|0), (1) 

where a, b are bosonic annihilation operators of the 
states. This does not have the appearance of possess¬ 
ing any entanglement. In first quantized notation, one 
would however write 

{ipa{ri)il;b{r2) + i>a{r2)ipb{ri)), ( 2 ) 
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Understanding entanglement in indistinguishable sys¬ 
tems is important to understand not only from a fun¬ 
damental perspective, but increasingly in the context 
of technological applications. The resource capabilities 
in quantum optical systems has been a topic of long¬ 
standing investigation, due to their highly coherent prop¬ 
erties making them a natural candidate for quantum in¬ 
formation tasks H- Recent experimental advances in 
coherent control of Bose-Einstein condensates (BECs) 
have opened the door for utilizing these systems for tasks 
such as quantum metrology [2^, interferometry 2^, 
quantum simulation and quantum computing 2M- 
In particular, the ability of realizing BECs on atom chips 
allows for the possibility of excellent coherent control , 
and performing tasks such as spin squeezing [^ . This 
allows for a naturally scalable architecture where multi¬ 
ple BECs could be produced on the same chip, allowing 
for the possibility of entanglement between BECs [IMi. 
While entanglement between two BEC has not been real¬ 
ized to date, entanglement between a BEC and an atom 
has been demonstrated [3l|. While entanglement crite¬ 
ria for optical systems has been relatively well-studied 
[3^ [3^ , the analogous case for BECs has not been per¬ 
formed to the same level Q. 


which has the appearance of an entangled state. The 
amount of entanglement in such systems that is opera¬ 
tionally accessible is difficult to define and has been a sub¬ 
ject of a long debate. Initially the debate was motivated 
by the problem of nonclassical properties of a two-mode 
single photon field Later many propositions for 

defining entanglement of indistinguishable particles were 
proposed [5l4l3l|. which were concluded by showing, that 
only a part of entanglement comin g fr om symmetrization 
can be operationally extracted [lj-[l^ . Recently more ef¬ 
ficient way of extracting this kind of entan glem ent was 
proposed by a method called mode splitting [l9|. 


In this paper we consider a hybrid scenario, in which 
two systems of bosons constitute two physically distin¬ 
guishable components. For simplicity we assume that 
both components contain the same fixed number N of 
indistinguishable particles, which can be prepared in two 
different modes. In such a physical system, bipartite 
entanglement between the two components is well de¬ 
fined, and can be experimentally detected and used as 
a resource. The main goal of our paper is to present 
correlation-based criteria for detection of this kind of en¬ 
tanglement. A correlation based approach [ssl - lssj is de¬ 
sirable from an experimental point of view as they are 
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most closely related to measured observables. However 
it cannot be directly applied to the case of multi-boson 
states, because it involves measurements performed on 
single localized particles. In order to extend applicability 
of correlation-based entanglement detetection to multi¬ 
boson systems, we define the measurement basis for the 
two components, which enables us to treat the entire sys¬ 
tem as two effective (A^-l-l)-level systems. Having defined 
the basis we apply quadratic correlation tensor criterion 
from [ 3 ^, which allows for the efficient detection of bi¬ 
partite entanglement between the two components. Note 
that the above method of probing entanglement utilizes 
the entire structure of subsystems that can be defined 
within the components. Furthermore, we demonstrate 
that using this entire structure is a necessary condition 
for detecting entanglement - the compound observables 
with higher-dimensional eigenspaces (e.g. total angular 
momentum observables) lead to weaker and unsatisfac¬ 
tory conditions. Despite its simplicity, the model pre¬ 
sented in this paper can be used to characterize entan¬ 
glement in two experimentally significant physical sys- 
tems: entangled spinor Bose-Einstein condensates 
and non-classical light created during spontaneous para¬ 
metric down conversion (SPDC) [i^. We choose our ex¬ 
amples based on realistic methods that could be used 
to generate entanglement in BECs and optical systems, 
which provide a test-bed for examining macroscopic en¬ 
tanglement [ 3 ^. l4lj|. 

This paper is structured as follows. In Sec. [TTlwe intro¬ 
duce our approach for bipartite bosonic entanglement. In 
Sec. IHIl we give various examples of our scheme applied to 
specific states, including the maximally entangled state, 
parametrically down converted state, two BECs evolving 
under a interaction, and two BECs in a two-mode 

squeezed state. In Sec. 113 we contrast our approach to 
alternative methods based on spin operators. We finally 
show our conclusions in Sec. El 


II. ENTANGLEMENT CONDITION 

We hrst recall the simplest version of the geometric 
entanglement condition [35l| in the context of two d- 
dimensional systems. Arbitrary quantum state of the 
entire system can be uniquely described by its correla¬ 
tion tensor T, which is defined in the Gell-Mann matrix 
basis in a natural way [l^ 

d-1 

p = -^(^0 0 Mo+Y^ T^jMi 0 Mj ( 3 ) 

d^-i 

-b (d - I) ^ {To^Mo 0M,+ T,oMi ® Mo)), 
where the coefhcients are 

^ 2(d- ® 


where Mi {i = 0,..., — I) are the generalized Gell- 

Mann matrices (see e.g. [^) and Mo = , ^ - 1. 

V d{d—l)/2 

The Gell-Mann matrices gene ralize the previous ap¬ 
proach considered in Refs. which were more fo¬ 

cused on the SU(2) case to SU(d). 

The correlation tensors are elements of a real vector 
space with a usual scalar product 

{X,Y) = J2x,,Y,,. ( 5 ) 

id 

A state p is separable if it can be expressed as a convex 
sum of product states 

Psep = ^ ', PiPA ® PB^ (®) 

i 

with Pi > 0 and Pi = 1- In the language of the correla¬ 
tion tensors this decomposition reads 
where 0 and represents a single 

qudit state. The separable correlation tensors form a 
convex set. In a consequence, when p is endowed by T, 
one has the following implication: 

p is separable =b d^^prod (T, > (T, T), (7) 

or, equivalently, 

max(r, < (T,T) =b p is entangled. (8) 

'Y'pTod 

The proof of the above implications is presented in [35l| . 
The maximization of the left-hand side of m is given 
by the highest Schmidt coefficient, Tmax, of the tensor T 
|44l| . Therefore, the condition 

T 

e = —- > 1 =b p is entangled, (9) 

J- max 

where T = {T,T) = '^ijT^j, is a simple entanglement 
witness. Note that the determination of T„iax requires 
knowledge of the entire tensor T and in general can be 
a hard computational problem. Since Tmax < I, we can 
rewrite the condition ® as 

T > I =b p is entangled. (10) 

Note that the original form of the entanglement crite¬ 
rion presented above is formulated only for the situation 
in which a subsystem is associated with a single d-level 
particle, and cannot be directly applied to the case sub¬ 
systems consist of many indistinguishable particles. 

Now we adapt the entanglement criterion to the case of 
bosonic systems. Let us consider the most general state 
of a two-component system of N two-mode bosons 

N N 

1^) = EE^^'1^- ® 1 ^- ObIOb, (11) 

k^o 1^0 

where denotes the i-particle state in the spatial mode 
X and some bosonic mode a. The states within given 
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spatial modes are given in particle number representa¬ 
tion, whereas the tensor product symbol between the 
spatial modes indicates full distinguishability between 
them. For simplicity, we will use the abbreviated no¬ 
tation \i)Wj)\ = omit the tensor product 

symbol. The core of our construction is the following 
dehnition of a new basis: 


|fc) = |iV-fc,fc), (12) 

where the vector |fc) is defined as the {k-\- l)-th eigenvec¬ 
tor of the last (diagonal) generalized Gell-Mann matrix of 
dimension fV-|-l. The definition (TT^ maps the two mode 
system with fixed number N of particles into a {N -1-1)- 
level system, which can be understood as a single logical 
qudit, where d = N -\-1. In the new basis the state OB 
takes the form: 

N N 

\'^) =^^Ckl\k)A\l)B, (13) 

fc=0 1=0 

and is formally equivalent to a state of two qucKts, with 
d = A^-|-l. As entanglement criterion we take the con¬ 
ditions or cni) with the summation indices i,j = 
{l,...,iV2 + 21V}. 

A specific example of (fT51) is the maximally entangled 
state of two qudits (with d = -|- 1 > 2), which can be 

written 


|V^W) = ^^f;|fc)^|fc)B. (14) 

The entanglement identifier for this state can be evalu¬ 
ated to be 


r(^)=iV-b2 (15) 

and Tmax = 1 Q- Since = A^ -|- 2 > 1, we see that 
it is entangled for all A^ > 0, as expected. 


III. EXAMPLE APPLICATIONS 


mediated cavity based methods 1^ . [30l| and interac¬ 
tions with state-dependent forces [IB. The type of en¬ 
tanglement that is produced by such an interaction has 
an exotic structure with fractal properties (the “devil’s 
crevasse”) [13, 

To produce the state, two spin coherent states are ini¬ 
tially polarized in the direction, and have a 
interaction applied to them. For BECs a general spin 
coherent state is written [s^ 

|a,^)) = ^(aat+/36t)^|0,0) (16) 

= (17) 

We define the Schwinger boson operators as 


=a'^h + h'^a 
= —ia^b + ib^a 

S^ = a^a-b^b. (18) 


After the entangling operation we obtain 

^ p-iSfSJr J_ J_\\ 1 1 


IV’BEc) = e 

N 


1 


vw 


N 


E' 


V 2 ' V2// V2’ V2 

gi(Af-2fe)r g-i(7V-2fe)r 


V2 


V2 


= hi\N-k,k)\N-l,l), 

k,l—0 


\N-k,k) 

(19) 


where bki = (^) (^)e 2 fc)(Ar 21 )t^ interac¬ 

tion is periodic in time r = 7r/2 hence we consider only 
T = [0,7r/2). 

When the number of particles per BEG is A = 1, the 
system is exactly equivalent to two qubits. In this case, 
the non-vanishing elements of the correlation tensor are 
equal to = I and = sin 2t. We may then 

evaluate: 


= l + 2s\n^2T ( 20 ) 


Here we present a few examples of applications of the 
condition (l8|) to detect entanglement of various classes 
of states. Our choices of examples are motivated by ex¬ 
perimental realizability. We show one example related to 
photonic systems showing parametric down conversion, 
and two examples for the BEG case involving an 
entangling gate and two-mode squeezing induced by en¬ 
tanglement with coherent light. 

A. Two BECs entangled with an interaction 

For our hrst example, we examine the case where two 
BECs are entangled using an interaction. Several 
proposals for such an interaction exist, including photon 


and 

T^l = 1- (21) 

We therefore find that > 1 for all r ^ 0. The max¬ 
imal value of = 3 is achieved for r = 7r/4, at which 
point we obtain a Bell state. All the expected behavior 
for an entanglement detector in this case is seen. 

For larger A, the complexity of the entanglement de¬ 
tector increases. For example, for A = 2, we obtain 

( 22 ) 

and 

tSL = ^ + 3^vE2, (23) 
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where c = 53 + 48 cos 4r + 24 cos 8r + 3 cos 16r. We there¬ 
fore illustrate our results numerically for TV < 8, as shown 
in Fig. [TJ In Fig. [T]^a) we plot which exceeds 

unity for all r ^ 0 (see Fig. [1]). While this shows that 
entanglement is present for all times by virtue of (ITOl) . 
the amount that this exceeds unity decreases with N. In 
fact the amount of entanglement, as measurement by the 
von Neumann entropy 

5 =-Tr(pilog2Pi) (24) 


where pi is the density matrix with the partial trace over 
BEC 2, is known to increase with TV. This is natural 
as the Hilbert space of the system increases with TV, 
hence the amount of entanglement can be also expected 
to grow. The maximally entangled state (HI has an en¬ 
tanglement equal to S'max = log 2 (A^ + !)• 

To remedy this behavior one may instead look at the 
stricter criterion given according to (IH]). Fig. mb) shows 
the behavior of with the entangling time. We see 
that again this exceeds unity for all times indicating the 
presence of entanglement. In contrast to we see 

that the amount the bound is exceeded by grows with 
TV, and more resembles the behavior of the von Neumann 
entropy. We may therefore observe that is a more 
sensitive indicator of entanglement, although it has the 
additional overhead of the necessity of calculating Tmax- 
We note that there is a very simple criterion that can 
be used to detect entanglement, restricted to the pure 
states and any particle number N. The observation is 
that: a pure state j'l/') is entangled if and only if the 
lengths of its local Bloch vectors are less than 1 [47| . For 
the BEC state, two Bloch vectors are of the same length 
given by: 


dW||2 =5^ + ^5^4 cos2"2Zt, (25) 


where 

^ 2Cat(TV)(iV + l)nU(^-^ + l) 

4^nU(^+*) 

^ (iV + l)2Cat(N) 1 

0 qATjy ^ 


(26) 

(27) 


and Cat is the Catalan number. For the first few values 
of TV, ||a(^)|p is equal to: 

||a(il|p = cos^ 2r (28) 


7 ( 2 ) 


7(4) 


i2 I 3 4 3 4 

f = - + - cos^ 2 t + — cos^ 4 t 
16 4 16 


7(3) 112 


15 1 1 

= — + - cos® 2t + - cos® 4r + — cos® 6t 
12 8 4 24 


,2 47 35 o 35 o 

5 5 

+ — cos® 6r + —— cos® 8t 
64 512 


It is easy to see that the length of the Bloch vectors 
||^(Af )||2 _ 2 ^ only for r = {0,7r/2}. For other values of r. 





FIG. 1. (Color online). Various entanglement identifiers (a) 
(b) (c) the von Neumann entropy S for two BECs 

entangled via a interaction as a function of r. 


||q(^)||2 ^ 2, which proves the entanglement of the BEC 
state. 


For TV = 1, there exists such r = 7r/4 that the state 
is maximally entangled (||a(^-’|P = 0). For TV = 3,4 the 
minimal value of | | p increases. For TV > 4 the length 

of the Bloch vectors is a decreasing function of TV and for 
TV —>■ oo, ||a('^^|p —>■ 0, which means that the BEC state 
again becomes maximally entangled (in this limit). 
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B. Two BECs entangled via two mode squeezing 


We next consider an entanglement generating proce¬ 
dure where we adapt experimental techniques used to 
entangle two atomic ensembles to BECs. Specifically, 
we consider the protocol experimentally realized by Juls- 
gaard et al. [i^ and proposed in Refs, In this 

protocol, one starts with two BECs, both polarized in 
the direction, which are illuminated by coherent light 
in a superposition of two polarization modes 


’^ 2 ’ 72 ^^’ 72 ’ 


(29) 


where the coherent state of light of amplitude a is 


1^) ^g-|a|V2ga/\/2(cWdt)|o)_ 


(30) 


and c, d denote the two polarization modes. The light 
and ensembles then interact via the ac Stark shift, as 
the light is detuned off the excited state resonance. The 
interaction Hamiltonian is 


= (^r + 5|) (31) 

where Sf ^ are the total spin operators as before, and 
J* = c^c — d^d is the spin operator for the light. This 
couples the three systems and produces entanglement of 
the type BEC-light-BEC. One then measures the light in 
the basis. This yields the state [5l| 

\t, Tic, rid) OC E ^ricTid {ki,k2)\ki,k2), (32) 

ki ,^2 

where, up to an irrelevant phase, 

An,,nAki,k 2 ) = sin”'= Ccos"-'7 

(33) 

and ^ = 2t{ki + k 2 — N) + j. The proportionality on 
the state (1321) refers to the fact that the state is unnor¬ 
malized as it is the state after a projective measurement 
has been made on the photons. The ric^d are the ran¬ 
dom measurement outcomes for the number of photons 
detected in the basis. 

We have performed numerical calculations for N up 
to 40 atoms. We assume parameters such that IctP = 
N, which is the typical experimental regime [^. As 
the photon number in a coherent state should obey ric + 
rid ~ |ap, we take ric = , rid = N - ric- Fig. [2] 

shows how the entanglement parameter depends on 
the interaction time t for various N. We find that there 
is a strong odd/even effect where the periodicity is 7r/4 
for even N and 7r/2 for odd N. This follows from the 
exponents in (1331) where for even N our choice gives ric = 
rid, and we have the simplification 

X cos"° (4t(fci -I- ^2 — A^)), (34) 




FIG. 2. (Color online). The values of for the macroscopic 
entanglement state (|32|l as a function of t for N odd and even, 
separately. 


which has an extra factor of two in the cosine. For odd 
N, we have an additional factor of cos^ in ([M|) which 
then has twice the periodicity of the even case. In ei¬ 
ther case we find that our approach successfully detects 
entanglement, with > 1 for all t. 

In the experimental system, in fact the typical 
timescale of entanglement is the short-time limit, with 
typical values t ~ 2.5/N [H, [5l|. We plot in Fig. 
[21 as a function of N. We see that the amount that the 
bound is exceeded does not generally increase with N, 
and it stays mostly in the vicinity of ~ 2. This is 
expected as in the short-time limit with interaction times 
in the region of t ^ 1/N the amount of entanglement is 
generally of the order of = 1 case and remains a con¬ 
stant. This behavior was observed also for the in¬ 

teractions where for interaction times of r ~ 1/7V the 
von Neumann entropy is of order unity. This originates 
from the dimensions of the ac Stark shift Hamiltonian 
(l3T|l which is of order 0{H) ~ N^, and therefore to keep 
the total evolution equivalent the time needs to be 

scaled as t ^ 1/N. 


C. State from parametric down conversion 

Our final example is the resulting multiphoton state 
in the spontaneous parametric down conversion process 
creating an entangled bright squeezed vacuum state ■ 
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FIG. 3. (Color online). The values of Tmax and 

for the macroscopic entanglement state (|32l) with t = 2.5/N. 
For each case, > 1 what proves entanglement. 


The state is 

I OO 

I^Z-PDc) =- KVN+1 (35) 

cosh K 

where 

1 ^ 

-m,m)A\N-m,rn)B, (36) 

^ m=0 


and |A^ — m,m)x denotes a state of N photons {N — m 
horizontally and m vertically polarized photons) in the 
spatial mode X {X = {A, B}). 

As the number of photons in the state (l35ll is not fixed, 
the condition (l8|) cannot be directly used to analyze en¬ 
tanglement properties. This can be circumvented by as¬ 
suming that we have photon number resolving detectors 
and therefore access to the number of detected photons. 
This allows us to post-select events that come from the 
maximally entangled state of an arbitrary number of pho¬ 
tons. After measurement, we have the state 

1 ^ 

|m)^|m)B, (37) 

^ m—0 

with probability 

Pn = -^^^tanh^"-A:(n-h 1). (38) 

cosh K 


For the states |'0^x) the value of condition ([5]) is equal 
eW = N + 2 (see Eq. (11^1) ). The average value of e over 
different number of emitted photons can be evaluated to 
be 


OO 

Savg = X] = 

N=0 


1 

16 cosh'^ K 


(39) 


X (15 cosh 2 AT-I- 6 cosh 4Ar-I- cosh 6 AT — 6), 


which is greater than 1 for all AT >0. We thus again find 
that the modified criterion can detect entanglement. We 
note that the quantity Eavg has only a statistical meaning 
and cannot be measured in a single run of the experiment. 


IV. GELL-MANN VERSUS SPIN OPERATOR 
BASIS 

In the criterion ([8]) we employed a decomposition to 
the space spanned by the Gell-Mann matrices (HI). One 
may consider an alternative approach based only on spin 
operators, which could be equally used to defining an 
entanglement criterion. We explain in this section the 
shortcomings of such an approach and the advantages of 
the general case ([8]). 

Let us define 

Tlj = ]^Tr(p5* 0 Sn, (40) 

where are the Schwinger boson operators (1151) . The 
condition for entanglement has the same form as m , ex¬ 
cept that now 

j-'iN) = Y. t'^ (41) 

i, 0 ={x,y,z} 

and T,^ax refers to the correlation tensor T'. While this is 
a perfectly valid entanglement criterion, we can examine 
why such an approach may not be as useful as the general 
case As an example, let us consider the BEG state 
(US discussed above. Its correlation tensor T' has the 
following non-vanishing elements 



nx = cos^^-\2r), 

(42) 


T'yz = T'zy = cos^"^(2r) sin(2r). 

(43) 

This yields 




= - cos^^“^ 2 t(cos^'^ 2r -|- sin^ At). 

(44) 


After diagonalization the correlation tensor T' has the 
form —T^y}, and T^ax is simply equal to 

max(r,^,j., \TyJ). We plot the results of using such a cri¬ 
terion in Fig. 0] for various N. While for N = 1 such a 
criterion works (i.e e > 1) for all r, for N = 2,3 the 
criterion only works in restricted regions of t. For N > 3 
it completely fails and one cannot reveal entanglement 
for all r. This is in contrast to the Gell-Mann approach 
(Fig. [T]) where entanglement can be detected for all N. 

The reason why we observe the difference between the 
approaches utilizing the Gell-Mann matrices and the spin 
operators is that three Schwinger operators {Sx,Sy and 
Sz) do not span the entire space of the dimension d = 
V-|-1. Moreover, if we express the Schwinger operators in 
the number state basis, we observe that each component 
of the spin-^ operator can be expressed in terms of N{= 
d — 1) Gell-Mann matrices: 

- d—1 

~ 2 ^ ~ -^|(2d(fc-i)-fc2+fc+2)5 (45) 

k^l 
^ d—1 

^ 2 ^ ~ M d(d~i) ^ fc(fc-i) ^ (46) 

k^l 

S'=5 (47) 
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FIG. 4. (Color online). The values of e parameter for 
N < 20. If e > 1, entanglement of the BEC state can be 
revealed by means of the condition based on the spin opera¬ 
tors. 


j=^(M4 + M5), (48) 

-{M-j + v^Afg). 

This means that in order to describe all three compo¬ 
nents of the spin operator S we need only 3(c? — 1) out 
of — I Gell-Mann operators. However, considering 
the higher powers of the s pin operators one can build 
a complete operator basis and use it to express all 
Gell-Mann matrices (see Appendix Thus while the 
criterion based on spin operators can in principle detect 
entanglement, higher powers of the spin operators would 
be in general necessary to detect entanglement. Using 
only the linear terms of the spin operators appears to 
miss out in capturing the necessary ingredients to show 
entanglement is present, at least for the cases that were 
examined in this study. 


For example, for d = 3 


= 


sy = 


= 



V. CONCLUSIONS 

We constructed a correlation-based criterion for de¬ 
tecting entanglement between two separate (distinguish¬ 
able) systems of bosons, which are indistinguishable 
within each subsystem. This is achieved by a change of 
basis from the standard occupation number basis to the 
effective N + 1 dimensional one, which can be treated as 
the eigenbasis of the N + 1-dimensional diagonal Gell- 
Mann matrix. This operation allows for the first time 


the application of the entire theory of entanglement de¬ 
tection based on correlation tensors. We have applied 
our theory to several realistic experimental examples, 
such as entangled BEGs, a light-induced two-mode 
squeezed BEGs, and photons from parametric down con¬ 
version. We find that in each case our criterion effectively 
detects entanglement for the full range of parameters. 

In order to evaluate the criterion in complete general¬ 
ity, a full set of Gell-Mann correlations would be required, 
which would be a challenging task experimentally. It is 
possible to derive an analogous criterion m based only 
on spin operators which would considerably simplify the 
detection process. Unfortunately this criterion appears 
to be highly insensitive, failing to detect entanglement 
altogether for large particle systems. This would suggest 
that correlations beyond lowest order in the spin oper¬ 
ators are necessary to detect entanglement. Using the 
full set of Gell-Mann correlations gives a much stronger 
signal of entanglement, and thus this approach appears a 
more robust method of detecting entanglement. One ad¬ 
vantage that our scheme has over other approaches (e.g. 
using negativity 0) is that it is based on correlations, 
rather than tomographically reconstructing the full den¬ 
sity matrix. Thus once a sufficient number of correlations 
are measured, it gives a direct proof of entanglement, 
even if the complete set is not measured. Exactly what 
correlations out of the full Gell-Mann matrix set capture 
a given type of entanglement is a more difficult task and 
is an open question. 
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Appendix A: The Gell-Mann matrices for d = 3 in 
terms of components of the spin operator 


Let S^,Sy and be the components of the spin-1 
operator presented in Eq. (1481) . The set of the three- 
dimensional Gell-Mann matrices can be expressed in 







terms of those operators in the following way: 


Ml 

M 2 

M 3 

Mi 


0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

i 

0 


1 0^ 
0 0 
0 0 y 

0 l^ 
0 0 
0 0 y 

0 0^ 
0 1 
1 oy 

-i 0 
0 0 
0 0 


V2 

{s^f - {syf, 

^2 

= -^{sy + 8^8^ + s^sy), 
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